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ABSTRACT 

Neutron star magnetic fields may have polar caps (PC) that are offset from the 
dipole axis, through field-line sweepback near the light cylinder or non-symmetric 
currents within the star. The effects of such offsets on electron-positron pair 
cascades are investigated, using simple models of dipole magnetic fields with 
small distortions that shift the PCs by different amounts or directions. Using a 
Monte Carlo pair cascade simulation, we explore the changes in the pair spectrum, 
multiplicity and energy flux across the PC, as well as the trends in pair flux and 
pair energy flux with spin-down luminosity, L^^- We also give an estimate of 
the distribution of heating flux from returning positrons on the PC for different 
offsets. We find that even modest offsets can produce significant increases in pair 
multiplicity, especially for pulsars that are near or beyond the pair death lines 
for centered PCs, primarily because of higher accelerating fields. Pair spectra 
cover several decades in energy, with the spectral range of millisecond pulsars 
(MSPs) two orders of magnitude higher than for normal pulsars, and PC offsets 
allow significant extension of all spectra to lower pair energies. We find that 
the total PC pair luminosity Lpair is proportional to Lgd, with Lpair ~ 
for normal pulsars and Lp^h ~ lO^^Lsd for MSPs. Remarkably, the total PC 
heating luminosity for even large offsets increases by less than a factor of two, 
even though the PC area increases by much larger factors, because most of the 
heating occurs near the magnetic axis. 



1. Introduction 

Electron-positron pair production by rotation-powered pulsars is believed to be essen- 
tial to the generation of charged particles in their magnetospheres and winds, as well as 
to the generation of their coherent pulsed radio emission. Pairs can be produced through 
electromagnetic cascades initiated by accelerated particles and several sites of acceleration 
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have been proposed. Near the neutron star surface above the magnetic polar caps (PCs), 
electrons may be accelerated either in vacuum gaps (Ruderman & Sutherland 1975) or in 
space-charge limited flow gaps (Arons & Scharlemann 1979). The gamma-ray photons radi- 
ated by the electrons will convert to electron-positron pairs by one-photon pair production 
in the strong magnetic fleld near the surface, and the pairs radiate synchrotron photons that 
produce further generations of pairs (Daugherty & Harding 1982). In outer gap vacuum 
accelerators, extending to near the light cylinder, seed photons can produce pairs through 
interaction with thermal X-ray photons from the hot neutron star surface to create a first 
generation of particles that accelerate through the gap radiating gamma-ray photons that 
produce further pairs (Cheng, Ho & Ruderman 1986). To date, all of the simulations of pul- 
sar pair cascades have assumed pure dipolc magnetic fields both to derive the accelerating 
electric field parallel to the magnetic field and to calculate the spectrum and multiplicity of 
the pairs. 

The configuration of neutron star (NS) magnetic fields is currently an unsolved problem. 
The original solution for a rotating magnetic dipole in vacuum (Deutsch 1955) shows that 
retardation causes the field lines near and outside the light cylinder to sweep back in a 
direction opposite to that of rotation, transitioning to an electromagnetic wave at large 
distances. The magnetic field sweepback causes the boundary of the open field volume, and 
the footpoints of field lines on this boundary that define the PC, to shift relative to the 
magnetic axis opposite to the rotation direction (Arendt & Eilek 1998, Dyks & Harding 
2004). The amount of the shift depends on the magnetic field inchnation to the rotation 
axis, with the largest offsets occurring for oblique rotators. In the last decade, numerical 
solutions of force-free magnetosphere models assuming ideal MHD (E ■ B = 0) (Contopoulos 
et al. 1999) show an even larger sweepback of the magnetic field near the light cylinder 
as it transitions to an MHD wind (Spitkovsky 2006, Timokhin 2006), as currents cause 
significant distortions of the magnetic field. The open field boundary and the PC of a force- 
free magnetosphere thus have a larger offset than the vacuum magnetosphere. Since the 
vacuum magnetosphere can accelerate particles but is not loaded with charges or currents 
and the ideal-MHD magnetosphere has charges and currents but cannot accelerate particles, 
a real pulsar magnetosphere lies between these two solutions. 

In addition to the guaranteed distortions of the dipole field caused by retardation and 
currents, it is possible that there are also distortions due to multipolar components near the 
neutron star surface. The emission from millisecond pulsars in particular shows evidence 
of non-dipolar field structure. The thermal X-ray pulse profiles of some millisecond pulsars 

(MSPs) show asymmetries which have been successfully modeled by including offsets of the 
emitting hot spot on the neutron star surface. Since the emission likely originates from PC 
heating, it is argued that MSPs such as PSR J0437-4715 (Bognanov et al. 2007) and PSR 
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J0030+0451 (Bogdanov & Grindlay 2009) have either offset dipoles or offset PCs. Modehng 
of the X-ray pulsations of the NSs in Low-Mass X-ray Binaries, the hkely progenitors of 
miUisecond pulsars, show possible evidence of even more extreme magnetic field distortions 
(Lamb et al, 2009), that could result from distortion of the global magnetic field by e.g. the 
crustal plate tectonics (Ruderman 1991). 

Wc have recently shown (Harding & Muslimov 2011, HMll) that distortions of a mag- 
netic dipole field that produce offset PCs lead to an increase in the multiplicity of PC pair 
cascades and to significant shifts of the death line for pair production of curvature radiation 
photons. We found that even offsets that are small fractions of a neutron star radius can 
enable high-multiplicity pair cascades in older pulsars that with pure dipole fields would be 
pair starved. There were a number of previous studies of how multipolar or non-dipolar neu- 
tron star fields could increase pair production (Ruderman & Sutherland 1975; Arons 1983, 
1997; Asseo & Khechinashvih 2002), but all of these focussed on the decrease in field line 
radius of curvature as the primary effect. Arons (1997), considering the decrease in field line 
radius of curvature provided by an offset dipole, found that offsets that are large fractions of 
a stellar radius, (0.7 - 0.8)i?ns; are required to move the pair death line enough to include the 
entire radio pulsar population. Medin & Lai (2010) proposed an even larger offset of 0.95 
i?ns, again only considering the decreased field line radius of curvature. Kantor & Tsygan 
(2003, 2004) considered the space-charge limited flow (SCLF) solution of a dipole field 
with the addition of a displaced dipole of smaller magnitude, that effectively produces nar- 
row, bent fiux tubes. They take into account both the increase in E'y and decrease in radius 
of curvature to derive lower death lines for pair production, and find that displacements of 
0.1it!ns for the smaller dipole component can significantly change the death lines for both 
curvature radiation (CR) and inverse Compton scattering (ICS) pair production. Our study 
also takes into account the increase in E\\ and particle acceleration energy, which is by far 
the strongest effect on pair multiplicity of offset dipoles or PCs. 

The study of HMll was based on derivation of the electric field parallel to the magnetic 
field near the PC assuming there is free fiow of particles from the neutron star surface, the 
SCLF flow boundary condition (Arons & Scharlemann 1979). The resulting current den- 
sity of primary charges in this case is equal to J = p^jC where p^j is the Goldreich- Julian 
charge density. This current distribution differs from the current density distribution de- 
manded by the global ideal-MHD models, so it is presently not clear whether the steady 
SCLF acceleration model will be compatible with a global pulsar magnetosphere. However, 
more realistic global magnetosphere models with resistivity and dissipation are in devel- 
opment (Li, Spitkovsky & Tchekhovskoy 2011). The current density distributions of such 
pulsar magnetospheres with dissipation are signiflcantly different from those with ideal-MHD 
(Kalapotharakos et al. 2011), with the ability to possibly adjust to the J ~ p^^c of the PC 
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cascades they rely on for their charge supply. Furthermore, Timokhin & Arons (2011) have 
performed particle-in-cell simulations of PC pair cascades and show that SCLF acceleration 
is stable for currents within a 10% range of p^^c. 

In this paper we extend the study of HMll to include a more general expression for the 
distorted magnetic field, treating two cases in which the caps from opposite magnetic poles 
are offset in the same (symmetric) or opposite (asymmetric) directions. In Section 2, we 
give the expressions for the magnetic fields in these two well as the corresponding 

Goldreich- Julian charge density, field line equations and parallel electric fields. Section 3 
describes the pair cascade simulation used to derive the results for pair multiplicity and fiux 
distributions over the PC, and pair spectra given in Section 4. We also examine the changes 
in the distribution of PC heating luminosity with offset PCs as well as the total pair and 
heating luminosity. 



2. Magnetic Field Configuration 

Derivation of a general non-dipolar magnetic field of a neutron star using a multipole 
expansion gives quite complex expressions (e.g. Assco & Khechinashvili 2002). Because 
of the mathematical complexity, this kind of model representation of the magnetic field 
is impractical. Since our main purpose in this paper is to derive an expression for the 
accelerating electric field near the neutron star surface, we require field expressions from 
which we can derive analytic formulae for charge density for input to Poisson's Equation for 
the electric potential. Wc therefore propose the following heuristic model of a non-dipolar 
magnetic field that provides distortions leading to offset PCs. We introduce an azimuthal 
asymmetry to the field lines of an originally symmetric dipole such that the field lines over 
half of the PC have relatively smaller radius of curvature and over the other half of the 
PC have larger radius of curvature. This can be done by adding an azimuthal dependence, 
controlled by the parameter a — £sin0, to the magnetic colatitude of Bj. and Bq, as in 
equations (1) and (2). Then, from the solenoidality condition, V • -B = 0, we can derive the 
azimuthal component of the distorted magnetic field, which is proportional to the parameter 
e. The resulting exact expression for the azimuthal component is rather cumbersome and can 
be simplified by expanding the appropriate terms in e. In our study, for the sake of simplicity, 
we retained only the lowest-order (linear) terms in e. Note that equations (1) and (2) can 
be used to achieve a significant distortion of the field lines while keeping the dipole axis 
inside the offset polar cap boundary, which significantly facilitates the subsequent analytical 
treatment. Thus, equations (1) and (2) represent global open field lines of a centered dipole 
that are pushed and bent toward the dipole axis on one side and pushed and bent away from 
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the dipole axis on the other side of the polar cap, respectively (see Figure 1). Consequently, 
one side of the PC is larger and the PC is effectively shifted from the center of symmetry. 

Case A: (the magnetic field is symmetric w.r.t. the magnetic equator) 

In this case, in magnetic spherical polar coordinates (r^, 9, (p) the magnetic field reads 



B ^ — 



r cos 6 + - {1 + a) sin 6 — cj) e sin 6* cos 6 sin(0 — 0o) 



(1) 



Case B: (the magnetic field is asymmetric w.r.t. the magnetic equator) 
In this case the magnetic field can be presented as 



r cos[6'(l + a)] + ]-0 sin[6'(l + a)] - ^ (j) e {e + smO cos 9) sin(0 - 0o) 



(2) 



where Bq is the surface magnetic field strength at the magnetic pole, r] = r/R^g is the 
dimensionless radial coordinate in units of stellar radius, i?ns, a = e cos(0 — 0o) is the 
parameter characterizing the distortion of polar field lines, and 0o is the magnetic azimuthal 
angle defining the meridional plane of the offset PC. We must point out that B as presented 
by formula (|2]) is solenoidal only approximately, since in the i?(^-component we have neglected 
corrections of order and higher than ~ {sOY ^^"^ -^r- and i?e-components we have neglected 
the corrections of order and higher than ~ (sO)^ and ~ (f^^)^, respectively. Figure 1 shows 
the field lines projected in the x-z plane for cases A and B. 

The corresponding Goldreich- Julian charge density is 

Case A: 



Pgj 



n Bo 



COS X + 2 



1 - 



COS( 



£— sm0sm( 



sin X sin 6 cos 6 



(3) 



Case B: 



PGJ 



Q Bo 
27rcT]'^ 



1 - ^ ) cos X + - 



K 



- ( 1 H ^ ) sin sin 



1 cos 9 cos 0+ 

2r]'^ 

sin X sin 9 



(4) 



where Q is the rotation rate of the NS, x is the pulsar obliquity, k ^ 0.15 I45/RI (-^45 — 
J/IO^^ g-cm^, i?6 = -Rns/10^ cm, / is NS moment of inertia) is the parameter accounting for 
the general- relativistic frame dragging (Muslimov & Tsygan 1992). 
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2.1. Electric Field in the Polar Cap Region 

If the magnetic axis lies in the x-z plane, a = e cos to account for an effective offset of the 
dipole axis in the x-z (00 = 0) plane and a = esincj) for an effective offset of the dipole axis 
in the y-z plane (00 = 7r/2), where < e < 1. By using formulae ([1]) and ([2]) we can get the 
equations of the magnetic field line, 

Case A: 

d = sin-^ , (5) 

Case B: 

e = ^ sin-i (/ e , (6) 

respectively, where x = r/Ri^c is the radial distance in units of the light cylinder radius, 
i?Lc = c/fi;0<^<lis the colatitude of a footpoint of polar field line normalized by 
the colatitude of the PC boundary. Also, in formula (Q f{a) is a parameter (^1) that can 
be more reliably determined via numerical tracing of the last open field line. For a rough 
estimate one can derive the analytic expression. 



/ 



sm 



TT 



2(1 + a) 



i(l+a) 



(7) 



that is applicable for a > and can be used to plot the field lines with x, y > 0. The field 
lines with x, y < can be plotted by using the coordinate reflection, x{y) — )■ —x{—y) and 
z — )■ —z. 

The fleld line radii of curvature (in units of Rlc) are 
Case A: 



Xr = —, r- a; 2 



(i-a) [l-(3 + a)(l-a)A/4f/^ 



3(1 + a)^ l + a/3- [1 -2a(l + a/2)/3]A/2' ^' 

Case B: 

" 3/e l + a/3-/2A/2' 

respectively, where 

A = x^+\ (10) 



Neglecting the static general-relativistic corrections, we can write the accelerating elec- 
tric fleld, assuming the boundary conditions of SCLF (cf. Harding & Muslimov 1998; here- 
after HM98). For 7] - 1 > 6*0, where Oq = sin"^[((]i?ns/c)^/^] is the standard PC angle. 
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Case A: 



3 ^ (1 -O COSX+ e2A ^(1 -^sinx 



Case B: 



Eu 



1 ^ / 



2 \l + a 



K 



3 e. 



3 ^ eiB (1 - ^ ) cos X + -— ^(1 - sin X 



8^ 



where Eq = {VtR^s/ c)'^{Br / B)Bq, and 



111 



(12) 



eiA = eiB = 1 + - a{rf - 1), 



e2A = a <! (1 + OCOS0 + /3 ^^^^^^^^ 



2a + 



5 — 3a 



7] 



{5-a)/2 



and 



where 



e2B = a 



/ 



1 + a 



[1 + COS(/)+- /3 £ 



a = X 



a/2 



1 + 3a 1 



3^(l+a)/2 



and [3 



[I + cos (/)o - - ^ cos(2(/) - 0o) 
5 



(13) 
(14) 

(15) 
(16) 



Now we shall present the expressions for E\\ that are applicable at the small altitudes, 
z = ri-l<dQ (cf. formula [19] of Harding k Muslimov 2001; hereafter HMOl), 



^ —3 £q (^9, 



a(l) 




4k cos X Gi{ku z) ^^[^'P. + 



i=l 

oo 



^ 9q {XOq j sin X cos 



i=l 



(17) 



where Gj = 1 — exp(— 7^2;) and Gj = 1 — exp{—^iz); A = 1 and A = /(1)/[1 + a(l)] in case A 
and B, respectively; /i = 1 + a(l) + [5 — a(l)]fi;/2; Jq and Ji are the Bessel functions, and fcj 
and ki are the positive zeros of Jq and Ji, respectively; 7 = ki/ XO^'^"'^^^ and 7 = fcj/A^^Q^"*-^''; 
and a(l) = e cos [0(1) — (po], where 0(1) = (f){ri = 1). The expression for the parallel electric 
field in Eqn ( [T71) above corresponds to the expression for a pure dipole magnetic field given 
in Eqn (19) of HMOl. The main differences here are inclusion of the azimuthal asymmetry 
incorporated in the a parameter, and the neglect of the static GR corrections. So if one 
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takes a = 0, Eqn will be equivalent to Eqn (19) in HMOl with H{1) = 5(1) = /(I) = 1 
and e = 0. 

Note that in formula flTTl) (p ~ 0(1) ^"^^ ^''^^ related via the field-line equation, 

" "^'^ = iT.tlli^al-X] ~ "^'^ ^ ezsin^[0(l) - M- (18) 

We see that the magnetic field configuration of Eqns (1) and (2) produces an E\\ that is 
significantly larger on the offset side of the PC, by the factor 6q'^"'^^\ than the E\\ for the case 
of a pure dipole field. This increase in E\\ results from the strong dependence on the field-line 
curvature of the Goldreich- Julian charge density in the solution of Poisson's equation for the 
electric potential (see also Kantor & Tsygan 2003). This dependence is embodied in the 
factor 6q in the pure dipole case, and in the factor Ol""^^^ in the offset PC case, so that the 
PC boundary is dependent on azimuthal angle. For the case of an offset PC, the effective 
PC angle, ^q"''^'*, is larger on one side of the PC, producing a larger E^, and smaller on the 
opposite side, producing a smaller Eh. 



3. Pair Cascade Simulation 

We model the electromagnetic pair cascades above a pulsar PC using a hybrid steady- 
state acceleration/Monte Carlo pair cascade simulation. The first stage of the calculation, 
based on the self-limited acceleration model of HM98, follows a primary electron that starts 
at the neutron star surface at magnetic colatitude ^ and azimuth with Lorentz factor 
7 = 1. An initial guess is made for the height he of the pair formation front (PFF) above 
the surface and this distance is divided into 500 equal linear steps of length As. The PFF is 
the height where the first pair is produced and we assume that the £"11 is completely screened 
above this point. This is a good assumption for several reasons. First, the electric field arises 
due to a small imbalance between the actual charge density and the local, rot at ion- induced, 
Goldreich- Julian charge density. It therefore does not require much additional charge to 
short-out this field. Second, the onset of pair cascading occurs very quickly (Daugherty & 
Harding 1996), so that the number of pairs produced per primary particle increases rapidly 
over small distances. Thus, as found in Arons (1983), the width of the PFF (the screening 
distance of the electric field) is very small compared to other dimensions of the problem. 
The electron advances in steps of As, gaining energy through electrostatic acceleration and 
losing energy through radiative losses. The PFF results from pairs produced by 7-rays of 
energy emin (in units of mc^) radiated by particles of energy 7min- The height of the PFF is 
then: 

Sc = min[S'a(7min) + Sp{e„,i,,)] (19) 
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where S'a (7min) is the distance required to accelerate the particle until it can radiate a photon 
of energy emin, and S'p(emin) is the pair attenuation length of the photon. The acceleration 
distance, Sai'jmm), is determined by first integrating the equation of motion of the particle 
to determine its energy as a function of its pathlength s: 

where E\\ is the electric field induced parallel to the magnetic field, 'j^^ and are the loss 
rates for ICS and CR. Expressions for and j^^ and details of their derivation are given in 
HM98. The pair production attenuation length of photons radiated by the particle through 
either ICS or CR is then determined. The attenuation length Sp{e), defined to be the path 
length over which the optical depth is unity, is given by 

r(e) = / " Tpp(ekB, e)ds = 1 , (21) 



where ds is the pathlength differential along the photon momentum vector k , Tpp is the 
attenuation coefficient for one-photon pair production and 6'kB is the angle between k and 
the local magnetic field direction. 

The method we use to compute the electron-positron pair attenuation length of the 
photons, Sp{e), has been described in detail in Harding, Baring & Gonthier (1997). Using 
equation fl^ . Sp{e) is computed by integrating the pair production attenuation coefficient 
of the photon along its path through the dipole field. The photon is assumed to pair produce 
at the point where r(e) = 1, and Sp{e) is then set to that path length. The two main inputs 
needed are the energy of the photon and its angle to the magnetic field, 6'kB- The energies of 
the radiated photons for the purpose of computing the PFF are taken to be (following HM98) 
ep = 13ecR for B < 0.1 B^r, and = 4ecR/3 for B > 0.1 B^r, where ecR = (3/2)(A/pc)7^ is 
the critical CR energy, A = A/27r is the electron classical radius, Bcr = 4.413 x 10^'^ G, and 
Pc = XcRlc is the field line radius of curvature (Eqn [8] and [9]). For the angle of the radiated 
photons at the emission point, we assume that 6*^3 = 2/7 for Compton scattered photons 
(Dermer 1990) and ^kB = for CR photons. The CR and ICS photons are assumed to be 
completely polarized in the parallel mode. Due to the curvature of the field lines, ^kB will 
grow as the photon propagates, roughly as sin6'kB ~ s/pc. To evaluate the pair production 
attenuation coefficient at each point along the path of the photon, we Lorentz-transform 
the photon energy and local magnetic field to the frame in which the photon propagates 
perpendicular to the local field. This is the center-of momentum frame for the created 
pair, where the attenuation has its simplest form and the photon energy is e^.,^ = esin6'kB- 
The one-photon pair attenuation coefficient is considered in two regimes. For B < 0.1i?cr, 
photon pair produce far above threshold, where the asymptotic expression in the limit of 
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large numbers of kinematically available pair Landau states (Tsai & Erber 1974, Daugherty 
& Harding 1983 [DH83]) can be used: 

= i^'A||,±(x), (22) 

f (0.31,0.15) exp(-|-) x«l 
H±ix) ^ < (23) 
[ (0.72, 0.48) x»l 

where x = /^cm/2-B', a is the fine-structure constant, B' = B/Bct is the dipole field 
strength at point s along the photon path. When B > 0.1 B^r, pair production will occur 
near threshold, where the above expression is not accurate. We thus include the factor 
/ = 1 + 0.42e~^'^ in x, introduced by DH83, as an approximation to the near-threshold 
attenuation coefficient. In this paper we compute the attenuation length averaged over 
photon polarization. 

The path of each input photon is traced through the magnetic field, accumulating the 
survival probability for pair production, Psurv^ 

Psurv(s) =exp|-r(s)} (24) 

Each photon may pair produce or escape, based on a combination of the running survival 
probability for pair production. For this simulation, we have neglected photon splitting since 
pair production dominates the photon attenuation in magnetic fields B < 10^^ G (Baring & 
Harding 2001). 

The "first guess" value of he, and thus also of Sc, sets the initial acceleration length, 
equation (120|) is integrated in discrete steps upward from the starting point, computing E\\ 
from equation f lT7|) . •y^,-. and j^^^ at each step. At each step, the pair attenuation lengths, 
Sp{e), of CR test photons radiated by the particle of energy 7(3) are computed from equation 
( 12T|) . The pair attenuation length, and thus the value of Sc, also computed at every step, 
is initially infinite, because the energy of the photons is small, but decreases with distance 
as the energy of the radiated photons increases. Although the photon attenuation length 
continues to decrease, the particle acceleration length is increasing and Sc has a minimum. 
This minimum value of he = {Sc)mm is adopted as the new height of the electron PFF. The 
electron is accelerated again with the new value of he, producing a new PFF at the next 
value of he = (5'c)min- The process is repeated, converging to a self-consistent value of he- 

Once the PFF height is established, the full pair cascade spectrum is computed using 
the Monte Carlo simulation, adopting the electron Lorentz factor as a function of height. 
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7(s), from the self-consistent PFF calculation for the acceleration region from s = to 
s = he- Above the PFF, it is assumed that £"11 = due to pair screening (HMOl), so 
that the electrons lose energy every step due to CR and ICS. The step size for s > he is 
limited to less than the distance over which they would lose 10% of their energy. Each 
step, the electron radiates a CR spectrum, divided into a number of logarithmically spaced 
bins. The number of CR photons in each energy bin, ncR, is determined by the energy 
loss rate and average energy in that bin. A representative photon from each bin having the 
average bin energy is followed to its pair conversion or escape point, as described above. 
The pair produced by the photon, or the escaping photon number, is then weighted by ticr- 
If a photon pair produces, the total energy. Landau state and parallel momentum of the 
electron and positron arc determined. Each member of the pair is assumed to have half the 
energy and the same direction of the parent photon. Each member of the pair occupying 
an excited state emits a sequence of cyclotron or synchrotron photons. The method used 
to simulate the cyclotron/synchrotron emission is similar to that of Daugherty & Harding 
(1996). If the particle Landau level is larger than 20, the high-energy limit of the quantum 
synchrotron transition rate (Sokolov & Ternov 1968) is used, in which case we assume that 
the photons are emitted perpendicular to the magnetic field in the particle rest frame (high- 
energy hmit). When the particle Landau level is smaller than 20, the exact QED cyclotron 
transition rate (Harding & Prcccc 1987) is used, in which case the angles of the emitted 
photons are sampled from a distribution. In both cases, the emitted photon polarizations 
are sampled from the corresponding polarization distributions. Each emitted photon is 
propagated through the magnetic field from its emission point until it pair produces or 
escapes. The cyclotron/synchrotron emission sequence continues until each particle reaches 
the ground state. By use of a recursive routine that is called upon the emission of each 
photon, we can follow an arbitrary number of pair generations. The cascade continues until 
all photons from each branch have escaped. The cascade pairs are binned in energy and the 
magnetic colatitude ^ and azimuth (p on the PC of the initial test electron. 



4. Results 

4.1. Pair Multiplicity and Flux Distributions 

For a given set of input parameters that inchidc pulsar period P, surface magnetic field 
Bq, inclination angle x, and offset degree e, direction 0o and symmetry (case A or B), the 
cascade calculation produces the pair spectrum per primary electron as a function of ^ and 
surface 0, dN^{E,^,(f))/dE, where we will express the pair energy E in units of mc^. The 
pair multiplicity, the number of pairs per primary electron, as a function of ^ and (p across 
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the PC is 

M^i,,,).r-'-W^,E (26) 

^ ^min 

while the total pair multiplicity from each PC is 

MPC=/ dc^ M^{i,cl>)el^Rl,idi. (26) 
Jo Jo 

Here, the PC angle varies with 0, 

sin^PC = (i?nsl^/c)(^+'^)/^ (27) 

which comes from setting ^ = 1 and r = i?ns in Eqn dS]). The total flux of pairs emerging 
from each PC is 



iV+(e,0)= / 



where 

ciiv+(E,e,0) rfiv+(E,e, 



iVp(e,0) (29) 



dE dE 
and where Np{^, 0) is the primary flux over the PC, 

iVp(e, 0) = ncj C7ri?2^ sin^ ^pc. (30) 

where ncj = Pgj{v = l)/e is the Goldreich- Julian number density at the NS surface. In Eqn 
( 15U]1 . we use the Goldreich- Julian charge density given in Eqns and dl]). Similarly, the 
energy flux of pairs from each PC is 

^,«,,).r-;tof)^,^ (31) 

^ ^min 

We have computed dN^{E,C,,(f))/dE for a range of parameters to explore the distri- 
butions of M_|_(^,0), iV+(^,0) and E^{C,,(f)) across the PC. Our results are dependent on 
neutron star equation of state (EOS) through the E\\ dependence on the frame-dragging, 
speciflcally through the k parameter, and on the neutron star radius, through the PC angle. 
There have been a number of measurements of neutron star mass for MSPs that give values 
well above the canonical IAMq and as high as 2Mq (Demorest et al. 2010). We have there- 
fore assumed different EoS for normal (non-recycled) pulsars and MSPs. For non-recycled 
pulsars, we take Mns = 1.45 and radius i?ns = 10 km, which give a moment of inertia 
/ = 1.13 X lO^^gcm^ (Lattimer & Prakash 2007). We take a rotating NS model (Friedman 
et al. 1986) for MSPs with M„s = 2.15 Mq, R^, = 9.9 km and / = 1.56 x lO^^gcm^. 
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Examples of the pair multiplicity and pair flux distributions over the PC are shown in 
Figures 2 and 3, for a pulsar with P = 0.3 s and Bq = 3 x 10^^ G, and for Case B with 
(j)Q = Tc/2 for different values of the offset degree e. A pulsar with these parameter values 
lies above but near the CR pair death line for e = 0, in which case the cascade produces 
a modest pair multiplicity. As shown in Figure 2, the maximum pair multiplicity M+(^,0) 
occurs for ^ ~ 0.5 and is nearly symmetric in 0, with only a small asymmetry caused by 
the (f) dependence of E\\. As e increases, the peak pair multiplicity increases but only in a 
small region of the PC toward the offset, in this case at = 270°, while the multiplicity 
decreases in the other side of the PC in the = 90° direction. For larger offsets, the region 
of high multiplicity expands as the PC boundary expands in the direction of the offset. The 
distribution of pair flux for this same case, shown in Figure 3, is somewhat different from 
the distribution of pair multiplicity. For e = 0, E^{C,,(f)) peaks in a small ring around the 
center of the PC and decreases outside the ring, again with a small asymmetry due to the 
dependence of E\\. The concentration of i?+(^,0) closer to the PC center is due to the fact 
that the Lorentz factor of the primary electron is maximum near the magnetic axis because 
the decrease in radius of curvature of the field lines raises the altitude of the PFF (see Figures 
4 and 5 of HM98). For e > 0, the pair energy flux increases and becomes asymmetric, but 
with a greater concentration toward the PC center than the multiplicity. 

Figure 4 shows the distribution of pair multiplicity M+{^, 0) for the case of a millisecond 
pulsar with P = 2 ms, i?o = 5 x 10^ G, and for Case A with 0o = 7r/2 for e between and 
0.6. These parameters place the pulsar below the CR pair death line with very small pair 
multiplicity for e = 0, although the distribution of M_|_(,^,0) can be seen to show a more 
pronounced asymmetry toward the 'favorably curved' field lines at cos > than in the case 
of the longer period pulsar shown in Figure 2. For pulsars with shorter periods, the second 
term in E\\ depending on 6o sin x cos is relatively larger since the value of is larger. As e 
increases, the peak pair multiplicity again increases in a small region of the PC toward the 
offset, but with a stronger asymmetry toward the cos0 > side of the PC. For millisecond 
periods there is also less increase of ellipticity of the PC with increasing e. The period 
dependence of the PC ellipticity can easily be seen from theh expression Eqn fl271) for the 
PC angle. If we define the elongation factor e = 6'pc(0 = 3tt/2)/6q = 6^^, then e is higher 
for smaller ^o- 

Figure 5 illustrates the effect of inclination angle on the distribution of pair multiplicity 
and pair flux for a millisecond pulsar with P = 3 ms, i?o = 4 x 10^ G, and for Case A with 
00 = tt/2 and e = 0.4. M+(^, 0) peaks in roughly the same region of the PC for x = 30° and 
X = 60°, but the asymmetry toward favorably curved field lines is is stronger in the case of 
X = 60°. Both the maximum multiplicity and flux are higher for x = 60°. 
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The values of peak multiplicity, M+(^ = 0.5, = 270°), for the radio pulsar population 
in the P-P plane are shown in the contour plots in Figures 6-8. Figure 6 shows contours of 
log(M_|_) in the case of a centered PC (e = 0), where we have assumed the different neutron 
star EoS for normal (non-recycled) and MSPs described above. It can be seen that M+ is high 
for short periods and higher P, reaching a maximum of a few times 10"^, and drops sharply 
towards the death line even on a log scale. About half of the normal pulsar population and 
most of the MSP population lies below the death hne, a result that has been noted in a 
number of previous papers (e.g. Arons 1997, Zhang, Harding & Mushmov 2000, Hibschman 
& Arons 2001, Harding & Mushmov 2002). In Figure 7 and 8, we show contours of log(M+) 
for offset PCs with e — 0.4 having symmetric (case A) and asymmetric (case B) offsets. 
In both symmetric and asymmetric offsets, the regions of high pair multiplicity spread to 
lower P and longer periods. The increases in are most dramatic for that part of the 
population near and below the pair death line having very low M+ in the centered PC case. 
The pair death line thus moves down through nearly the entire population of both recycled 
and non-recycled pulsars for an offset oi e = 0.4. For pulsars with high P, M_|_ changes much 
less with even a large offset, with multiplicity saturating below ~ 5 x 10^. This saturation 
of M+, noted by HMll, is caused by several effects. Even in the centered PC case, when 
the magnetic field increases above ~ 5 x 10^^ G pairs are increasingly produced in the low- 
lying Landau states (Baring & Harding 2001), which results in fewer synchrotron photons 
and fewer generations in the cascade. For increasing PC offsets, the resulting increase in 
accelerating field and particle energy initially produces higher M+ pair cascades, but as the 
increasingly energetic particles produce higher energy CR photons, the pairs are produced 
at smaller angles to the magnetic field. This results in higher average pair energies with 
pairs in lower Landau states, further reduces the number of synchrotron photons. There are 
some significant differences in the log(M+) contours for cases A and B. The death line in the 
normal pulsar population is lower for case B and overall values of M+ are a bit higher. For 
MSPs, the case B death line moves up for very short periods. 

The death lines for both £ = and e — 0.4 are flatter for the shorter period MSPs 
and the death hues at all periods are flatter as for £ = 0.4. As noted by HMll, particle 
acceleration in pulsars that require a large fraction of the open field voltage to produce pairs 
is limited by CR reaction. In this regime, the electron Lorentz factors to not reach the full 
voltage drop between the neutron star and the PFF, but instead reach a lower steady-state 
Lorentz factor as the acceleration is balanced by the CR loss rate. Even for centered PCs, 
acceleration in MSPs operates completely in the CR reaction limit (Luo et al. 2000, Harding 
& al. 2002). 
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4.2. Pair Spectra 

The total spectrum of cascade pairs, integrated over the whole PC, is 

Figure 9 shows the integrated PC pair spectra for fig = 3 x 12 G and a range of periods 
typical of non-recycled pulsars, for centered and offset PCs. The spectra exhibit turnovers 
at low and high energies that depend on period. For centered PCs, the low-energy turnovers 
occur around Lorentz factor E = 100 and the spectra extend up to ii^ ~ 10^, with the range 
increasing for shorter P. For offset PCs, the spectra extend to lower energies, by about a 
decade for e — 0.4. Thus the total PC pair flux increases because pairs can be produced 
at lower energies. Pair spectra for parameters typical of MSPs are shown in Figure 10. For 
both centered and offset PCs, the MSP pair energies are much higher that those of normal 
pulsars, by a factor of about 100. The spectra extend from a low-energy turnover at ~ 10^ 
to a high energy cutoff around lO'' for e = 0, and from ~ 2 — 3 x 10^ to ~ 10^ for e = 0.6. 
Thus the highest-energy pairs are nearly as energetic as the primary electrons! The large 
difference between normal and MSP pair spectra is due to the difference in field strengths. 
In lower fields, photons must pair produce at higher energies, decreasing the relative M_|_, 
even though the larger PC sizes of MSPs give smaller curvature radii. But as is the case for 
normal pulsars, the increase in for offset PCs results from lower possible pair energies. 



4.3. Total Polar Cap Pair Flux and Luminosity 

The pair flux emerging from each PC is 

(33) 

*^ -^min 

which we examine as a function of pulsar spin-down luminosity, Lsd- Figure 11 plots the PC 
pair flux vs. Lsd for the non-recycled pulsar population for different degrees of offset. The 
lines are least-squares fits to the points for each case, resulting in the following approximate 
expressions: 

e = (34) 
e = 0.2 (35) 
e = 0.4 (36) 

where Lgd.ss = Lgd/lO^^ ergs""^. The pair flux is closely proportional to spin-down luminosity 
for £ = 0, with the dependence somewhat flattening for offset PCs due to a saturation at 



N ■ 

^ ' pair 



1.6 X W^s-^Ll-^l^, 
5.5 X 10^4 s-iLO^Ifg, 
3.5 X 10^5 s-^L0^3^5, 
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high Lsd- Figure 12 shows the equivalent plot for MSPs with the least-squares fits giving the 
following relations between iVpair and Lg^' 

7Vpair= 9.5 X 10=^^3-1^0^5^^ (37) 

8.5 X 10=^^ s-^ L'^Jl^, e = 0.2 (38) 
3.1 X 10=^^s-i LO^^, e = 0.6 (39) 

The total pair luminosity from each PC is 

(40) 

Plots of the PC pair luminosity vs. Lgd for non-recycled pulsars are shown in Figure 13, with 



the least-squared fits yielding the relations: 

Lp,i,= 2.0x10^1 erg 8-1^0^365, e = Q (41) 

2.7 X lO^iergs-^LO^Og, £ = 0.2 (42) 

8.2 X 10^1 ergs-i L^^^^, e = 0.4. (43) 



The pair luminosity seems to also be roughly proportional to spin-down luminosity, with 
the dependence flattening for higher s, with the efficiency for pair luminosity being repair = 

Vir/^sd ~ 2 - 8 X 10-1 

The Lpair vs. Lsd for MSPs is shown in Figure 14, with the corresponding relations from 
the least-squares fits: 

Lpai,= 3.1 X 10=^2 g^gg-l ^0^6^^ (44) 

3.2 X 10=^2 ej.gs-1 LOd%, e = 0.2 (45) 
7.8 X 10=^2 ergs-i L^^^, s = 0.6. (46) 

The pair luminosity efficiency for MSPs is much higher than for non-recycled pulsars, with 
Vpair ~ 3 — 8 X 10"^, about an order of magnitude higher. Since MSPs have much lower 
surface magnetic field strengths, they need to use a higher fraction of the total open-field 
line voltage to produce pairs, resulting in a higher Typair- MSPs also have higher 7-ray 
luminosity efficiencies (Abdo et al. 2009), because the primary particles which radiate the 
highest energy emission reach higher energies before pairs screen the accelerating field. 



4.4. Polar Cap Heating 



Pair cascades from particle accelerators near the PC will inevitably produce backflowing 
particles in the process of pair screening of the electric field. The backflowing particles will 
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accelerate through the same potential drop as the primary particles and deposit this energy 
on the neutron star surface, increasing the surface temperature of the PC. The radiation 
from the hot PC is predicted to be in the X-ray band (Arons 1981, HMOl), and observations 
have provided evidence for emission from hot PCs, especially in the case of middle-aged 
and millisecond pulsars (Zavlin 2007). The predicted X-ray luminosity from heated PCs 
in the case of SCLF accelerators with centered PCs (HMOl) roughly agrees with observed 
luminosities. With the possibility of offset PCs increasing the pair multiplicity near the 
neutron star surface, it is important to check that the returning positron luminosity does 
not overheat the PC and violate the X-ray luminosity constraints. For this estimate, we 
do not perform a detailed screening calculation with pair dynamics to derive the screening 
scale length, as in HMOl. Rather we will compute the maximum returning positron fraction 
which will give a conservative limit on the positron heating luminosity. 

The maximum fraction of returning positrons from the PFF can be estimated as (see formula 
[33] in HMOl) 

P+(^pff) 1 



PGj(2:pFp) 2 



\ _ P(^pff) 

PGj(2;ppp) 



(47) 



where Zppp = r^ppp — 1 is the PFF dimensionless altitude and ^(^^ppp), the primary charge 
density at the PFF, is derived by inserting 77 = 1 inside the curly brackets in Eqns (|3]) and 
dl]) and r] = r/ppp outside the brackets. 

The heating flux over the PC from returning positrons can then be estimated as 



dS 



cp+(2;ppp,^,0)7(^PPF,'C,0) (4J 



where 7(-Zppp, ^, 0) is the primary electron energy at the PFF on the field line with surface ^ 
and 0. The reason that we use the the primary energy instead of the potential drop between 
the surface and the PFF to estimate the energy of the returning positrons at the surface is 
that for some cases the particles become radiation-reaction limited before reaching the PFF, 
so that the potential drop would be an overestimate. 

Figure 15 shows the distribution of positron heating flux over the PC, dL+{^, (f))/dS, for 
the case of a non-recycled pulsar with P = 0.1 s and Bq = 3 x 10^^ G, for different e values, 
while Figure 16 shows results for a MSP with P = 2 ms and B = 2 x 10^ G. It is evident 
that most of the heating flux is concentrated at the center of the PC near the magnetic axis 
near ^ = 0, the contrast between the PC center and edges being about a factor of ten. The 
returning positron fraction is maximum both near ^ = and ^ = 1, where the PFF is at 
the highest altitude. However, the '-/{Zp^^,^, (p) is maximum near ^ = 0, causing a maximum 
heating flux near the PC center. For e > 0, most of the offset side of the PC is heated with 
lower flux than near the PC center. There is also an asymmetry such that the favorably 
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curved field lines have higher heating flux, also because of higher returning positron energies. 
The substantially lower heating in the outer part of the offset side of the PC, which has a 
much larger area, keeps the increase in total PC heating due to the offset to a minimum. In 
fact in Figure 15 and 16, the physical size of the hottest area is comparable for all values 
of e that we examined. Thus the most intensely heated area is several times smaller than 
the canonical PC size, even for large e values. In the case of the MSP in Figure 16, this 
effectively smaller heated spot is shifted from the PC center by a smaller fraction of the 
canonical PC radius than the shift of the whole PC. For the non-recycled pulsar in Figure 15 
the heated area remains nearly centered while the whole PC is shifted by a larger fraction. 

The total heating power (of precipitating positrons) can be estimated as (cf. formula 
[61] of HMOl) 

L+^c p+{Zpp^,^,(f))-f{Zpp^,^,(f))dS, (49) 

where the integration is over the area of a sphere cut by the polar flux tube at the radial 
distance r^ppp. Figure 17 shows the dependence of L+, as a fraction of spin down luminosity, 
on pulsar characteristic age, r = P/2P, for different period and offsets. Comparing to the 
results in Fig. 7 of HMOl for the e = case, the L^/Lgd here is higher by factors of 2 - 
5 since the values we have given here are the upper limits on positron heating rate, while 
the HMOl results are from numerical computations of the i^y screening. Our estimates here 
nevertheless serve to test the heating of offset PCs against the observational constraints. 
In the case of both normal pulsars and MSPs, the increase in heating for e — 0.2 is very 
minimal. Even for e — 0.4, the heating luminosity increases by no more than a factor of 2, 
and for MSPs there is almost no increase except in cases where there are different degrees 
of pair screening for increasing offsets (e.g. P = 2 ms at large ages). 

Thermal components have been detected from both middle-aged pulsars and MSPs 
in the soft X-ray band. In the case of middle-aged pulsars like PSR B 1055-52 and PSR 
B06556-I-14, typically both hot and cool thermal components appear in the spectra. The 
hot components, with luminosity Lh, may be due to PC heating and the cool components 
from neutron star cooling, but in any case hot components provide a limit to any theoretical 
heating of the PC. For these pulsars, the PC heating efficiency L^/Lsd = 1-4 x 10^^ for 
B0656+14, and 5 x 10"^ for B1055-52 (DeLuca et al. 2005), which are within a factor of 2 of 
the our estimated maximum heating efficiencies. In the case of MSPs, several have measured 
thermal components, including PSR J0437-4715, J2124-3358 and J0030+0451. Since these 
pulsars are too old to have significant cooling, their thermal emission is very likely from PC 
heating. Their measured PC heating efficiencies, Lh/ Lgd ~ 6 x 10~^ (Zavlin 2007) are a 
factor of 2-3 lower than the maximum estimated efficiencies in Figure 17. 
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We have examined the change in magnetic flux over the PC as e increases. For an offset 
of £ = 0.2 for a 0.1 s pulsar, the magnetic flux on the offset side of the PC is typicafly a 
factor of 2 higher than for the centered PC, while on the side opposite the offset the flux is a 
factor of 3 lower, and for larger offsets the flux ratio can be several orders of magnitude. For 
MSPs, the contrast in flux is much smaller with ratios of only a factor of 4-5 in flux across 
the PC for offset as high as e = 0.6. 

5. Discussion 

We have investigated the effects of offset polar caps on electron acceleration and pair 
cascades near the neutron star surface, including the cases where the PCs of the two hemi- 
spheres are offset symmetrically by the same amount or asymmetrically. The asymmetric 
offset case would apply to the PCs shifted by retardation and/or currents of the global 
magnetosphere, while the symmetric case could apply to neutron stars with some interior 
current distortions that produce multipolar components near the surface. The asymmetric 
PC offsets have now been shown to be standard in pulsar magnetosphere global geometry. 
The PCs in the retarded vacuum magnetosphere are shifted in a direction opposite to that 
of the rotation by an amount that varies from a maximum of 20% of the standard PC radius 
at X = 90° to no shift for x = 0° (Dyks & Harding 2004). In terms of our offset parameter 
£, the ratio of offset to standard PC angle is 9^q/9q ~ 9^^ ^ so that for the vacuum case 
£vac — 0.03 — 0.1, where the large values of £vac apply to MSPs with large 9q. In the force-free 
ideal-MHD magnetosphere, the offset percentages range from 45% at % = 90° to 30% at 
X = 30° (Bai & Spitkovsky 2010), giving a range £imhd — 0.09 — 0.2, again with the larger 
values of £imhd applying to MSPs. 

From our calculations, we flnd that the pair multiplicity and pair flux for offset PCs is 
distributed very asymmetrically over the PC, with higher values on the side of the PC toward 
the offset. The regions of peak pair flux and multiplicity occur for £ = at about half the 
PC radius, symmetrically around the magnetic axis. As e increases, small regions of higher 
peak flux and multiplicity appear toward the offset and grow larger with increasing offset. 
The increase in peak multiplicity moves the pair death lines downward in P-P space, to 
encompass nearly the entire radio pulsar population for e = 0.4. Since such offset values are 
higher than those resulting from retardation and currents in pulsar magnetosphere models, 
they would likely require interior currents that produce large-scale non-dipolar flelds near the 
neutron star surface. Examination of the pair spectra reveal that the increase in multiplicity 
in pulsars with offset PCs comes from the extension of the spectra to lower pair energies. 

The discovery of pulsed gamma-ray emission from a large number of millisecond pulsars 
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by the Large Area Telescope (LAT) on the Fermi Gamma-Ray Space Telescope (Abdo et al. 
2009) has revealed light curves that are best modeled by narrow radiation gaps in the outer 
magnetosphere (Venter et al. 2009, Abdo ct al. 2010). Such narrow gaps require screening 
of the accelerating electric field over most of the magnetosphere by a pair multiplicity that 
is orders of magnitude larger than standard models of PC acceleration with no offsets are 
able to produce (Harding & Muslimov 2002). Since PC offsets that are only small fractions 
of a stellar radius result in large increases in pair multiplicity, and there is evidence for such 
offsets in MSPs, offset PCs may be a viable explanation for the larger-than-predicted pair 
activity in MSPs. 

We find that the integrated pair fiux and luminosity from each PC is roughly propor- 
tional to the pulsar spin down luminosity, both in the case of non-recycled pulsars and MSPs 
with MSPs being somewhat more efficient in converting spin-down to pair luminosity. For 
pulsars with offset PCs, both pair flux and luminosity increase with offset but to a greater 
degree for low spin down power. Due to a saturation of the pair multiplicity, pulsars with 
high Lsd have only modest increases in PC pair fiux even for large offsets. We estimate that 
the Crab pulsar produces a pair fiux from each PC of about 10^* pairs s~^ in the case of no 
offset, ~ 2 X 10^*^ pairs s^^ for e = 0.2 and ~ 5 x 10^^ pairs s^^ for e = 0.4. The fiux from 
both PCs is still more than an order of magnitude smaller than the pair fiux required to 
account for the radiation from the nebula, which is estimated to be about ~ 4 x 10^° pairs s~^ 
(DeJager et al. 1996). The trend of pair flux and pair luminosity proportional to spin down 
power may have interesting implications for explaining observed trends in pulsar pulsed and 
un-pulsed X-ray luminosity that are also proportional to spin down power (Vink et al. 2011). 

1/2 

Such a trend contrasts that of observed 7-ray luminosities which arc proportional to 
(Thompson et al. 1997, Abdo et al. 2010b). The 7-ray luminosity trend can be understood 
if the emission is produced by the primary PC current (proportional to lI^^) accelerated by 
a constant voltage of about 10^'^ V (Harding 1981, Arons 1996). The X-ray luminosity trend 
may be understood if both the pulsed emission from the magnetosphere and the unpulsed 
emission from the pulsar wind nebula (PWN) are produced by pairs. It is generally believed 
that pairs are producing the emission from PWNe, and a number of high-energy emission 
models argue that secondary electron-positron pairs produce the pulsed X rays high in the 
magnetosphere through synchrotron and/or ICS (Cheng et al. 1986, Romani 1996, Takata 
et al. 2007, Hirotani 2008, Harding et al. 2008). 

The asymmetric distribution of particle acceleration and pair multiplicity that results 
from offset PCs should produce asymmetries in observed pulsar emission. The structure 
and energetics of the proposed slot gaps (SG) that form between the boundary of the open 
magnetic field and the upward curving PFF (Arons & Scharlemann 1979), and can accelerate 
particles to high altitude (Muslimov & Harding 2004), could be strongly affected by a an 
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offset PC. The particle Lorentz factor 7 in the SG, which is expected to reach curvature 
radiation-reaction hmit such that 7 cx E^^^'^, will be larger on one side of the PC, producing 
CR emission power proportional to E\\ that is larger than for a dipole field. In older pulsars 
that do not produce enough pair multiplicity and screening to form SGs in centered PCs, 
SGs may form on only one side of an offset PC. An azimuthal asymmetry of both the 
radiation power and width of the SG would change both the 7-ray luminosity and the 
sharp 7-ray peaks that are due to caustics formed by radiation from trailing edge field lines 
(Dyks & Rudak 2003). Since the PC offsets that result from retardation and currents in 
pulsar magnetospheres occur toward the trailing side of the PC, the 7-ray peaks should be 
enhanced relative to the off-peak emission that is due to emission along the leading-edge 
field lines. The ratio of E\\, and thus 7-ray fiux, between trailing and leading edges of the 
SG is predicted to be ~ + ^)/(l ~ ^)) which can be larger than an order of magnitude 

for pulsars having short periods. SG model light curves assuming the emission asymmetry 
predicted in vacuum or non-ideal MHD magnetospheres fit the pulsar light curves measured 
by Fermi significantly better that those of symmetric SGs (Harding et al. 2011, DeCcsar ct 
al. 2011). The asymmetry in pair flux in offset PCs should produce asymmetries in pulsar 
radio emission if the radio flux is proportional somehow to the pair flux. In fact, evidence 
for such asymmetries in pulsar radio emission have been observed in the form of partial cone 
emission (Mitra & Rankin et al. 2011). In a future study, we will explore the effect of an 
offset PCs on the acceleration of particles in the SG at high altitudes and the change to the 
shapes of 7-ray light curves. 
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Case-A field lines, c = 0.2 




Case-B field lines, e = 0.2 




Fig. 1. — Field lines of distorted magnetic dipole having an offset polar cap in the x-z plane 
and offset parameter e = 0.2 for the symmetric (case A, top) and asymmetric (case B, 
bottom) polar cap offsets. 
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Fig. 2. — Contours of log of pair multiplicity as a function of radial distance from the 
magnetic axis (in units of neutron star radius) and magnetic azimuth for P = 0.3 s, Bq = 
3 X 10^^ G and x = 60°, for different values of offset parameter e. 
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Fig. 3. — Contours of log of pair energy flux (in units of ergs"^ cm^^) as a function of radial 
distance from the magnetic axis (in units of neutron star radius) and magnetic azimuth for 
P = 0.3 s, So = 3 X 10^^ G and x = 60°, for different values of offset parameter e. 
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Fig. 4. — Contours of log of pair multiplicity as a function of radial distance from the 
magnetic axis (in units of neutron star radius) and magnetic azimuth for P = 2 ms, Bq = 
5 X 10* G and x = 45°, for different values of offset parameter e. 
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Fig. 5. — Contours of log of pair multiplicity (top) and pair flux (bottom in units of cm~^) 
as a function of radial distance from the magnetic axis (in units of neutron star radius) and 
magnetic azimuth for P = 3 ms, Bq = A x 10^ G and x = 45°, for e = 0.4 and different 
values of inclination angle x- 
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LOG (P) 

Fig. 6. — Contours of log of peak pair multiplicity as a function pulsar period P and period 
derivative (P), for a centered PC and inclination angle x = 60° for normal pulsars and 
X = 45° for MSPs. The NS radius i?ns; niass (in Solar mass units) and moment 
of inertia = J/lO^^gcm^ refer to different NS equations of state used for normal and 
millisecond pulsars and are described in the text. Radio pulsars with measured P from the 
ATNF catalog (Manchester et al. 2005, ,http: // www. atnf . csiro . au / research / pulsar / psrcat ) 
are plotted as black dots. 
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Pair Multiplicity, 8 = 0.4, case A 
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Fig. 7. — Same as Figure 6 but for a symmetrically offset PC with e = 0.4 
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Pair Multiplicity, 8 = 0.4, case B 
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Fig. 8. — Same as Figure 6 but for a asymmetrically offset PC with e = 0.4 



- 33 - 




E (mc^) 

Fig. 9. — Spectra of pairs (pairs/(s mc^)) from one polar cap, for a surface magnetic field 
strength i?o = 3 x 10^^ G, different periods and two different degrees of offset e. 



- 34 - 




Fig. 10. — Spectra of pairs (pairs/(s mc^)) from one polar cap, for a surface magnetic field 
strength Bq = 2 x 10^ G, different periods and two different degrees of offset e. 
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Fig. 11. — Total pair flux (in units of pairs s ^) from each polar cap as a function of spin 
down luminosity Lsd for non-recycled pulsars, for different values of offset parameter e. 
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Fig. 12. — Total pair flux (in units of pairs s ^) from each polar cap as a function of spin 
down luminosity Lsd for millisecond pulsars, for different values of offset parameter e. 
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Fig. 13. — Total pair luminosity (in units of erg s ^) from each polar Ccip clS db function of 
spin down luminosity Lsd for non-recycled pulsars, for different values of offset parameter e. 
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Fig. 14. — Total pair luminosity (in units of erg s ^) from each polar Ccip clS db function of 
spin down luminosity Lsd for millisecond pulsars, for different values of offset parameter e. 
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Fig. 15. — Contours of log of positron heating flux (in units of ergs~^ cm~^) as a function 
of radial distance from the magnetic axis (in units of neutron star radius) and magnetic 
azimuth for P = 0.1 s, -Bq = 3 x 10^^ G and x = 60°, for different values of offset parameter 
e. 
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Fig. 16. — Contours of log of positron heating flux (in units of ergs~^ cm~^) as a function 
of radial distance from the magnetic axis (in units of neutron star radius) and magnetic 
azimuth for P = 2 ms, i?o = 2 x 10^ G and x = 45°, for different values of offset parameter 
e. 



-41 - 




10'° 

10^ 10^ 10^ 10^ 10^ 10^ 10^ 10 



T (years) 



Fig. 17. — Positron heating luminosity as a fraction of spin down luminosity Lsd vs. pulsar 
characteristic age r, for different pulsar periods and values of offset parameter e 



